Non Thermal Equilibrium States of Closed Bipartite Systems 
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We investigate a two-level system in resonant contact with a larger environment. The environment 
typically is in a canonical state with a given temperature initially. Depending on the precise spectral 
structure of the environment and the type of coupling between both systems, the smaller part may 
relax to a canonical state with the same temperature as the environment (i.e. thermal relaxation) or 
to some other quasi equilibrium state (non thermal relaxation). The type of the (quasi) equilibrium 
state can be related to the distribution of certain properties of the energy eigenvectors of the 
total system. We examine these distributions for several abstract and concrete (spin environment) 
Hamiltonian systems, the significant aspect of these distributions can be related to the relative 
strength of local and interaction parts of the Hamiltonian. 

PACS numbers: 05.30.-d, 03.65.Yz, 75.10.Jm 



I. INTRODUCTION 

In a composite but closed quantum system in which 
a smaller central system S is weakly coupled to a larger 
environment C, most of the (pure) states of the total sys- 
tem for a given energy (and possibly some additional con- 
straints) exhibit properties of thermal equilibrium states 
with respect to the smaller part jlj, i.e. there exists a 
so-called dominant region in Hilbert space in which the 
entropy of the central system is close to its maximum 
value under the given constraints. Therefore for most 
pure initial states of the total system, the state of the 
central system shows decoherence and some kind of ther- 
malization; it typically approaches a quasi-equilibrium 
canonical state with a temperature given by the spectral 
properties of the environment Q • 

If the environment initially is in a thermal state with 
a given temperature and consists of many bands or of a 
broad continuum of levels, the central system typically 
relaxes to a thermal state with the same temperature. 
This type of relaxation will be called thermal relaxation 
in the following. Here we investigate to what extent cer- 
tain structures of the total system influence the reached 
(quasi) equilibrium state. We will relate this equilibrium 
state to the distribution of the energy eigenvectors of the 
system, or rather certain important aspects of this dis- 
tribution. We will show that there is a close relation 
between the two, and how this affects the equilibrium 
state for different system structures. 

We particularly focus on a single spin- 1/2 particle 
coupled to an environment of spin-1/2 particles. Re- 
cently, the properties of spin systems of different struc- 
ture (rings, stars, and others) have been subject of exten- 
sive interest. A lot of work has been done on the question 
of entanglement [1, 0) IE IE 0JH , their relaxation be- 
havior has been addressed [9|, L!JJ] and various techniques 
were suggested to make any spin interact with any other 
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FIG. 1: A two-level system S in contact with an environ- 
ment C consisting of two highly degenerate levels k, k' with 
degeneracies g k , g k i . 



spin [ll|,llj. 

Here we extend our analysis form our previous pa- 
per [l3[ regarding the controllability of relaxation behav- 
ior within these spin systems 



II. CANONICAL AND NON CANONICAL 
RELAXATION 

Figure [1] shows a two-level system (TLS) in reso- 
nant (<5s = 5c = 5) contact with an environment consist- 
ing of two "energy bands" k, k' of degeneracies g/., g^, 
respectively (for simplicity we use g = gk, g' — gk' in the 
following, typically g' > g). The coupling is assumed to 
be weak, the total system is described by the Hamilto- 
nian 

H = Hs + He + Hi nt . 

A non equilibrium state |l)s<£>|</>fc)c is depicted (here \(f>k) 
denotes an arbitrary pure environmental state in band k) . 

If this state is taken as the initial state of a Schrodinger 
time evolution of the total system, a relaxation to an 
equilibrium situation is expected in which the time- 
averaged reduced state operator of S is given by [l[ 
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which can be interpreted as a canonical state operator 
with inverse temperature 
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FIG. 2: z-component of the Bloch vector of S for the initial 
state depicted in figure [T] The environment consists of spins, 
curves are for different types of intra-environmental coupling. 
Fluctuations are due to the finite size of the system. The 
black line indicates canonical equilibrium. 



For a finite environment and weak random coupling, the 
reduced state of the central system after relaxation still 
fluctuates around ([T]), see 0. 

If the central system relaxes to state (P) for any pure 
state \<pk) of the environment in band k, it will relax to 
the same state for an initial state in which the environ- 
ment is completely mixed within band fc, po — 

lfe denotes the projector onto band k of the environment. 

Assume now that the environment is given by a large 
number N of two level systems with a homogeneous Zee- 
man splitting, Hq — ^2iLi8&z . The environment ini- 
tially is taken to be in a thermal state of temperature /3c, 

p = |l)(l|®ie-^. 

If each band of the environment separately leads to a 
relaxation into a state ([T]), the equilibrium state of the 
central system is given by the canonical state with /3s = 
Pc for large N. For finite N, the population of the excited 
state of S after relaxation is given by 



yi n iV + ll + e 5 ^ 7V + 1 jv»i 1 + 

Because of this we call the relaxation from an initial 
state |1} S <g> \<j>k)c t° tne reduced equilibrium state ([T]) 
canonical or thermal throughout this text. 

However, not all environments lead to canonical relax- 
ation of the central system. In [l3j we examined sev- 
eral types of spin environments and showed that many of 
these lead to an equilibrium state that differs from |T]). 
We call these deviations non canonical or non thermal. 
Figure shows the relaxation behavior for different types 
of environments. Obviously, not all relax to the same 
quasi equilibrium state. 



III. ENERGY EIGENVECTOR 
DISTRIBUTIONS 

We now correlate deviations from canonical relaxation 
with the distribution of the energy eigenvectors of the 
total system. We only consider the situation depicted in 
figure [TJ since environments with more bands in a canon- 
ical state can simply be derived from these results. As 
long as the interaction is weak relative to the band split- 
ting, the total system can be reduced to the subspace 
consisting of the "cross states" 

{ | Os, environment in band k ), 

| Is, environment in band k)}. (2) 

In the following we will always refer only to this (g + g')- 
dimensional subspace. 

Instead of the temperature we consider the population 
inversion Tr{<7 z /5s} of the central system for mathemati- 
cal convenience. The inversion of the canonical state ([T]) 
is given by 

I" \ - 9 ~ 9 ' 
^)c,n- g + g r 

The energy eigenvectors can be written in the form 

\e)^a e \0,Xe)+Pe\l,V s }, (3) 

where \Xe) is a state in band k! , and \rj e ) a state in band k 
of the environment. In the following we will use e as a 
discrete index running from 1 to g + g' to number the 
energy eigenvectors within the subspace of Hilbert space 
spanned by the states ([2]) . 

We expand the initial state po = <8> -y in terms 

of these eigenvectors, p — Yl e s' Po,ee' \ £ )( £ '\- Averaging 
over all times and tracing out the environment yields the 
equilibrium state of the central system, 

PS = J2 Tr c {|£>(s|} 

€ 

= ^|/3 e | 2 (|a £ | 2 |0)(0| + |/3 £ | 2 |l)(l|). 

€ 

The respective inversion is 

<^> = Ei&i 2 (i&i 2 -i^i 2 )- 

e 

We notice that 

A e = |/3 £ | 2 -|a e | 2 (4) 

is the inversion of |e). Since |/3 £ | 2 + \a e \ 2 — 1 we can 
rewrite the time-averaged inversion of the central system 
completely in terms of the A £ 's, 
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FIG. 3: A E -distribution for random Hermitian matrices from 
the GUE. The solid line shows the (normalized) probability 
density ([6| , the histogram is calculated form the eigenvectors 
of 400 different random matrices, g = 91, g — 364, AAJ: = 
2/1425 « 0.0014, <<r z ) can = -0.6. 



The first sum can be shown to be equal to g — g' , see 
appendix [AJ If we rewrite the second sum in terms of 
the variance AA 2 of the A e -distribution, we finally get 



(<?z) = (<7,) can + 
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The average deviation from the canonical equilibrium 
state is thus mainly given by the distribution of the re- 
duced states of the energy eigenvectors, in particular by 
its variance. 

As long as the width of the distribution is finite, there 
is always a deviation from the canonical inversion and 
therefore also from the canonical temperature. For a fi- 
nite environment this is always the case. We will now 
consider several types of environments and interactions. 



their eigenvectors. The parameters used are g = 91 
and g' = 364. 

The probability density for the A £ 's for random Her- 
mitian matrices from the GUE can be calculated analyt- 
ically and is given by 



P(A) oc (1 — A) s ' _1 (l + A) 
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(6) 



This result is derived in appendix [B] The solid line in 
figure [3] shows the distribution. 

The mean value of A e for this distribution is given by 



as expected from the general result in appendix |A1 This 
is again the mean inversion of the equilibrium state for 
canonical relaxation. 

The variance of this distribution is given by 
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d 2 (d+l) 



with d = g + g' . Since this is finite, the reached steady 
state will typically deviate from the canonical equilibrium 
state. 

For large systems, i.e. d —> 00 (for constant ratio g/g'), 
both AAg and (a z ) — (cr z ) can vanish. So in the thermody- 
namic limit, the quasi equilibrium reached equals canon- 
ical equilibrium. 



B. Random interaction 



bands 



IV. EIGENVECTORS DISTRIBUTION FOR 
DIFFERENT HAMILTONIANS 

A. Random Hamiltonian 

At first we assume that the two relevant bands in 
the environment are exactly degenerate and the transi- 
tions within S and C are exactly in resonance. In this 
case, Hs+Hc oc 1 (within the relevant subspace spanned 
by the states ©), thus we only need to deal with H- m t- 

If the coupling iJj nt between system and environment is 
modeled by a random Hermitian matrix with a uniform 
Gaussian distribution w(H int ) = Cexp(— ATr{_ff; nt } 2 ) 
(i.e. taken from the GUE [lH), the state of the central 
system typically relaxes to the expected equilibrium state 
under Schrodinger dynamics for the total system P, 0, 
[l3j |. The state fluctuates in time, the amplitude of these 
fluctuations decreases with the size of the environment. 

Figure [3] shows the A £ -distribution for these random 
Hermitian matrices. The histogram was obtained by 
choosing a number of random matrices with the given 
probability distribution and calculating the inversion of 



We will now discuss a system that is not completely 
random, but has a random energy exchanging coupling 
between the central system and the environment. If the 
environmental bands are strictly degenerate, the Hamil- 
tonian matrix has the form 
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(7) 



where A — 6 s — 5 C is the detuning between system and 
environment. The upper left block corresponds to the 
ground state of the central system, its therefore of di- 
mension g' , the lower right block (of dimension g) cor- 
responds to the excited level. These blocks are purely 
diagonal. 

The off diagonal block V (a g x g' matrix) , correspond- 
ing to energy exchange (canonical) coupling, is chosen 
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FIG. 4: Random interaction, degenerate bands. The two bars 
represent delta-peaks of different strength. AA^ = 6/25 = 
0.24. 



randomly with normalized Gaussian distributions for the 
real and imaginary parts of the matrix elements. 

For A ^ 0, system and environment are off resonance 
and the relaxation to the canonical equilibrium state is 
prohibited by energy conservation. So just A = is 
considered. 

Figure 2] shows the A e -distribution for a matrix of 
type (O, again for g = 91 and g' = 364. Instead 
of being single peaked, the distribution here is drasti- 
cally different and consists of delta peaks at A e = —1 
and A £ = 0, respectively. For a given Hamiltonian, there 
are exactly g' ' — g eigenstates with A e — — 1 and 2g eigen- 
states with A e — (assuming g' > g). The mean value 
still is A £ = (g — g')/(g + g'), as expected. 

However, the variance of the distribution is obviously 
considerably bigger. For the given parameters AA £ = 
0.24 as opposed to AA £ « 0.0014 for the completely ran- 
dom Hamiltonian. The deviation of the quasi equilibrium 
state for a Hamiltonian of type (|7|) from the canonical one 
is exactly 



{<Jz) - {°z) c 



z / can 



which is nicely demonstrated in figure[5j Since this result 
is independent of the system size, going to large environ- 
ments will not change the relaxation behavior other than 
reducing the amplitude of the fluctuations. Even in the 
thermodynamic limit the canonical equilibrium state is 
never reached. 
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FIG. 5: Time evolution of the z-component of the Bloch vec- 
tor of S for a random interaction and degenerate environ- 
mental bands, as described in section IIV B II Note that the 
canonical equilibrium state for the given parameters would be 
at (a z ) = -0.6. 
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FIG. 6: A TLS in contact with an environment consisting of 
two "bands" with equidistant level spacing. Both bands have 
the same width Se, the lower band consists of g levels, the 
upper band of g' levels. 
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2. Non degenerate bands 

The situation, again, changes if we introduce a finite 
spacing between the levels within each environmental 
band. For simplicity we will only consider equidistant 
levels and equal bandwidth 6e for both bands here. The 
lowest (and highest) levels of each band in the environ- 
ment are in resonance with the central system. This sys- 
tem is depicted in figure [5] (5s — 6c is considered in the 
following). The Hamiltonian matrix in subspace @ is 
given by. 



By introducing a small level splitting in the environ- 
ment (small compared to the interaction strength), the 
peaks in figure[3]get broader, especially the one at A e = 
gets flatter considerably and is stretched towards nega- 
tive A e . The variance of the distribution gets smaller. 
Figure [7] shows the A £ -distribution for a relatively small 
level spacing. 

When the level splitting is increased, the distributions 
becomes single peaked, with the peak close to its average 
and of similar height as the peak of the complete random 
Hamiltonian, see figure [5) A long tail towards higher A e 
prevails, therefore the variance is still considerably larger. 
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FIG. 7: Random interaction, equidistant spacing, small level 
spacing. AA^ « 0.19. The solid line shows the distribution 
for the complete random Hamiltonian. 




FIG. 8: Random interaction, equidistant spacing, larger level 
spacing. The solid line shows the distribution for the complete 



random Hamiltonian. AA e 



0.022. 
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FIG. 9: Schematics of the spin-star (left) and spin-ring (right) 
configurations. Typically, the environment consists of a lot 
more than five spins. 
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FIG. 10: A E -distribution for the spin-star configuration. 
14 environmental spins, the 2nd and 3rd excited bands are 
considered, corresponding to g — 91, g' — 364. AA^ » 0.216. 



system-environment interaction is 



N 



Hint = 



EE 



,M~(S) 

ij u i 



C. Spin environments 

We now consider a gapped spin or TLS as in figure Q] 
coupled to an array of spins, all with a Zeeman splitting 
equal to the central spin (or almost equal). Due to energy 
conservation, the system can be reduced to the situation 
shown in figure [1] for each pair of environmental bands. 
If the environment initially is in a canonical state, we can 
simply sum up over all bands, as long as the interaction 
between the environmental spins is small. 

The A £ -distributions of several spin-environments have 
been discussed in [l^ | , so we will only discuss them briefly 
here. 



for N environmental spins 
It has been shown in 11 



that if the coefficients 7 
are chosen randomly, the initial state depicted in fig- 
ure [4] typically does not relax to the canonical equilib- 
rium state. The Hamiltonian matrix in this case has the 
form ([?]), however with small fluctuations on the diag- 
onal, and the interaction part of the Hamiltonian ma- 
trix is only sparsely populated. Nevertheless, the A e - 
distribution shows some similarity to the one described 
in section lTV B 2l for small level splitting. Figure [TU] shows 
the A £ -distribution for this system. 



2. Intra- environmental coupling 



1. Spin-star configuration 

Figure [S] (left) shows schematically a spin-star config- 
uration, i.e. a central spin coupled to an array of envi- 
ronmental spins without mutual interaction. A typical 
environment should of course consist of a lot more than 
5 spins. The most general Hamiltonian describing the 



If mutual coupling between the environmental spins 
is introduced, the situation changes. Figure [S] (right) 
schematically shows next neighbor coupling in the envi- 
ronment (spin-ring configuration), but other configura- 
tions are possible as well. As long as this coupling is 
weak, the system can still be considered band wise. 

The interaction typically leads to a level splitting 
within the bands which in turn leads to a A e -distribution 
similar to the one described in section HV B 21 with bigger 
level splitting. Figure [11] shows the A e -distribution for 
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FIG. 11: A e -distribution for a spin-ring configuration. 14 en- 
vironmental spins, the 2nd and 3rd excited bands are consid- 
ered, corresponding to g — 91, g — 364. AAJ; » 0.0295. 
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FIG. 12: A e -distribution for a spin-star configuration with 
inhomogeneous Zeeman splitting of the environmental spins. 
14 environmental spins, the 2nd and 3rd excited bands are 
considered, corresponding to g = 91, g' — 364. AAJ; » 0.0548. 



a spin-ring configuration. The shown distribution is for 
a a x ® a x next neighbor coupling. The distributions for 
different kinds of coupling, e.g. Heisenberg coupling, are 
similar. 



3. Inhomogeneous Zeeman splitting 

If the individual environmental spins each have a dif- 
ferent Zeeman splitting, the situation becomes similar 
to the one discussed in section IIVB 21 Figure [12] shows 
the A e -distribution when the Zeeman splittings of the en- 
vironmental spins are homogeneously distributed within 
a certain range. The distribution again shows a peak 
around the mean value A £ = (o"z) can , although broader 
than in the previous case. 



WIDTH OF THE DISTRIBUTION AND 
SPECTRAL WIDTH 



of the environmental spectrum from zero to the minimal 
variance of the A e -distribution. A similar transition can 
be observed for many different environmental spectra. 

In order to relate different types of spectra we split the 
Hamiltonian matrix (in the considered subspace) in its re- 
spective diagonal and off diagonal parts. The off diagonal 
part H g (V and of 0) describes the interaction be- 
tween central system and environment, while the diago- 
nal part -ffdiag describes the environmental spectra alone, 
if system and environment are in resonance. The diago- 
nal part is always taken to be traceless. The quantity we 
use to compare different spectra is the relative strength 
of the environmental part to the interaction part, 



\ 



Tr(g d 2 iag 
Tr(i? 2 ff ) 



For a completely random matrix (GUE) this relation is 
given on average by 



Vi 



R.GUE 



which only depends on g/g', not on the actual size of the 
system. 

Figure [TBI shows the variance of the A E -distribution for 
three different types of environmental spectra. In all 
three cases the environment consists of 14 spins, the 2nd 
and 3rd excited bands are considered, g = 91, g' = 364, 
as described in section HV CI 

The solid line corresponds to the spin-ring config- 
uration as described in section IIV C 21 The intra- 
environmental interaction is taken as a a x ® a x + a y (£> a v 
next neighbor coupling, the central system is randomly 
coupled to each environmental spin. The dashed line cor- 
responds to the configuration described in section HV C 3[ 
there is no mutual interaction between the spins in the 
environment, but their Zeeman splitting is inhomoge- 
neous. The central system is again coupled randomly to 
each environmental spin. The environmental spectrum 
for the dotted line is the same as for the dashed line. 
However, the interaction between the central spin and 
each environmental spin is modeled by a x ® a x + <j v ® a y . 
The corresponding average values for a random matrix 
from the GUE are AA 2 = 2/1425 « 0.0014 and Vr,gue » 
1.46, respectively. 

We notice that for each type of coupling and environ- 
mental spectrum there is a distinct minimum of AA 2 for 
similar values of Ur close to, but not exactly at the av- 
erage value for the GUE matrices. 



Figures [4] [Jj and [8] indicate that there is a continuous 
transition from a situation far from canonical to an al- 
most canonical relaxation, depending on the environmen- 
tal spectrum. What has been changed is the "strength'' 



VI. CONCLUSION 

We have characterized situations under which non- 
thermal states should result as quasi-equilibrium states. 
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FIG. 13: Variance of the A e -distribution over the relative 
strength of the environmental spectrum for different types 
of this spectrum. 



For a spin- 1/2 particle weakly coupled to a larger en- 
vironment, there is a close relation between the quasi- 
cquilibrium state of the small quantum system coupled 
to a larger environment and the distribution of certain 
properties of the energy eigenvectors of the total system. 
The equilibrium state is directly given by the width of 
this distribution. The spectral structure of the environ- 
ment and the exact form of the coupling has a strong 
influence on the eigenvector distribution. To show this 
we have considered both abstract system Hamiltonians as 
well as Hamiltonians for structured spin environments. 

Furthermore, there is a close relation between the 
width of the A 6 -distribution and the strengths of both the 
local (diagonal) and the interaction (off diagonal) part of 
the Hamiltonian. By changing certain parameters within 
each system, there is a distinct minimum of the A £ -width 
for a value of the relative strength that's close to the re- 
spective value for GUE matrices. This relative strength 
can thus give an indication whether for a given system 
relaxation to or close to a thermal state can be expected 
without calculating the full A e -distribution. The relative 
width gives an indication how to choose the system pa- 
rameters properly to achieve a certain type of equilibrium 
situation. This should be of help when designing a spin 
system for special ("non-thermal") relaxation behavior. 

We thank the Deutsche Forschungsgemeinschaft for fi- 
nancial support. 



APPENDIX A: TOWARDS EQ. © 

Here we show that ^ £ X e — g — g'. 
^Ae = ^Tr s {afTr c | £ )( £ |} 

e e 

= J2 Tt {( & " ® iC ) \ £ )( £ \} = Tr U^ s ® iC )}> 

e 

the last equality follows from the fact that the energy 
eigenvectors are a complete orthonormal basis in Hilbert 
space. If we now use the basis {|0, to') , |1, to)} to calcu- 
late the trace (to and to' denote the levels in band fc 



and fc', respectively), we see that the g kets \l,m) 
yield l's, while the g 1 kets |0, to') yield — l's and the total 
trace equals g — g 1 ■ 



APPENDIX B: DERIVATION OF EQ. © 

If we introduce the basis {\k:m)} (1 < to < g^) for the 
band k and the basis {\k':m)} for the band k', we can 
write the state of the total system as (g = <?/., g' = g^') 

a' a 

m—l m— 1 

The reduced state of the central system becomes 



Ps 



= Tr c P=Y1 



0m I 



I>l m | 2 |l)<l| 



We are interested in the distribution of the inversion of 
the energy eigenstates of certain Hamiltonians. Since the 
inversion is determined by the population of each level, 
we will derive the distribution for the population of the 
ground state, p = XLlV'Oml 2 - 

For simplicity we introduce a single index n to label 
the amplitudes instead of the double index Oto or Ito. 
n runs from 1 to d = g + g' (the reduced Hilbert space 

dimension). In this notation po = X)n=ilV'ra| 2 - 

We now split the amplitudes into real and imaginary 
part, tp n — Xk + ixk+i, where 1 < k < 2g' corresponds 
to 1 0> <0 1 and2 ff '+l < fc < 2d corresponds to |1}(1|, there- 
fore po = 2fe=i x k- The combined probability density of 
the first 2g' amplitudes Xk for eigenvectors of random 
matrices from the GUE is given by [3| 



r(d- 9 ') 



2ff' s d-g'-l 



fc=l 



The desired probability density for the population of the 
ground state po is given by 



PpiPo) = / d 2s 'x P a (xi, . . . ,X2 9 >) 5(po 

^ fc=i 



4 i, 



integrating over the unit sphere in 2(/-dimensional space. 
The integral yields 



P P (Po) 



r(d) 



r(d-g> 



-s'-i 



Transforming to A = 1 — 2po finally gives the probability 
density ©, 



P(X) 
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